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Abstract: The Pauli groups are ubiquitous in quantum information theory because 
of their usefulness in describing quantum states and operations and their readily 
understood symmetry properties. In addition, the most well-understood quantum error 
correcting codes—stabilizer codes—are built using Pauli operators. The eigenstates of these 
operators—stabilizer states—display a structure (e.g., mutual orthogonality relationships) 
that has made them useful in examples of multi-qubit non-locality and contextuality. Here, 
we apply the graph-theoretical contextuality formalism of Cabello, Severini and Winter 
to sets of stabilizer states, with particular attention to the effect of generalizing two-level 
qubit systems to odd prime d-level qudit systems. While state-independent contextuality 
using two-qubit states does not generalize to qudits, we show explicitly how state-dependent 
contextuality associated with a Bell inequality does generalize. Along the way we 
note various structural properties of stabilizer states, with respect to their orthogonality 
relationships, which may be of independent interest. 

Keywords: contextuality; stabilizer; non-locality 


1. Introduction 

The Pauli operators are ubiquitous in quantum information theory, typically used as an operator basis 
to decompose multi-particle states or circuits. There are many combinatorial and geometrical structures 
that arise in finite-dimensional quantum mechanics that are intimately related with Pauli operators, and 
these structures are known to have favourable properties for use in quantum information processing. 
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Examples include Mutually Unbiased Bases (MUBs [1,2]), Weyl-Heisenberg covariant Symmetric 
Informationally-Complete Positive Operator Value Measures (SIC-POVMs [3]) and spherical/unitary 
designs [4,5]. While the qubit (two-level system) versions of thePauli operators are the most familiar, 
the generalization to qudit (('/-level systems) Pauli operators is mathematically straightforward. What is 
not straightforward, however, is whether construction techniques, like those in the works listed above, 
will generalize and retain their desirable properties when we switch to general qudit systems. 

With regard to fault-tolerant universal quantum computation (UQC), Pauli operators arise naturally 
within the context of stabilizer error-correcting codes and so they are intimately related to this task 
too. On the other hand, the Gottesman-Knill theorem tells us we cannot see better-than-classical 
computational performance using quantum circuits restricted to (i) operating on Pauli eigenstates, (ii) 
using Pauli measurements; and (iii) using gates that inter-convert Pauli operators (Clifford gates). In this 
model of fault-tolerant computation, the use of qudits appears to offer an advantage over qubits in terms 
of the efficiency associated with a magic state distillation routine [6]. 

In this work we examine whether sets of stabilizer states ( i.e ., Pauli eigenstates) exhibit quantum 
contextuality. More precisely, we examine whether known instances of quantum contextuality for 
two-qubit states generalize to two-qudit states when d is an odd prime. The most directly comparable 
work is by Planat [7] although he restricts his search for contextuality to the subset of stabilizer states 
that can be represented as real vectors and allows non-power-of-prime dimensions. 

Because d = 2 is the only even prime dimension, we sometimes see striking dissimilarities in 
combinatorial or geometrical structures, depending on whether the prime Hilbert space dimension is 
d = 2 or d > 2. An example of this difference can be seen in the context of a discrete Wigner function 
(DWF) where the associated geometrical structure is simplified in (power-of-) odd prime dimensions 
[8,9]. We will see a similar distinction when we seek to find instances of state-independent contextuality; 
such instances are fairly easy to concoct using multiple qubits but seem impossible to find using multiple 
odd-prime dimensional qudits. 

Nonlocality can be understood as a special type of contextuality, and so Bell inequalities can be recast 
in the graph-based contextuality formalism. We derive such a decomposition for a family of two-qudit 
Bell inequalities and discuss the related orthogonality graphs, graph parameters and state-dependent 
noncontextuality inequalities. 

We begin by providing the necessary definitions and mathematical background along with a 
discussion of quantum contextuality. The subsequent results section is partitioned into three main 
sections. The first two subsections relate to state-independent and state-dependent contextuality 
respectively. The remaining subsection provides some additional considerations related to graph-based 
contextuality that should be borne in mind, and discusses how these apply to our current investigation. 

2. Mathematical Preliminaries 

In order to appreciate the results of Section 3 we will first need to provide the necessary mathematical 
background and definitions. We begin by discussing the generalized Pauli group, the associated stabilizer 
states and the unitary operations that inter-convert them. In the next subsection we introduce quantum 
contextuality and its graph-theoretical interpretation. This latter aspect of our work requires a dedicated 
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subsection that lists all the purely graph-theoretical relationships and results that are used elsewhere in 
the paper. 


2.1. Stabilizer States and Clifford Gates 

2.1.1. . Stabilizer States 


Throughout, we always assume the dimension, d, of a single qudit to be a prime number. Because 
of the mathematics involved, most of the subscripts, superscripts and arithmetic will use elements 
from Z,/—the set of integers modulo d. Occasionally we will also need to use Z * d —the set of non-zero 
integers modulo d. Define the generalized Pauli “shift” and “phase” operators, respectively, as 

X\j) = \j + 1), Z\j) — u J \j) where a; = exp(2vr i/d) (1) 


so that in the d — 2 case these reduce to the familiar <j x and a z operators. Taking all products and powers 
of these X and Z operators produces the set of Pauli operators 


P{x\z) 


i xz X x Z z (d = 2) 
X X Z Z (d > 2) 


( 2 ) 


where we have introduced a phase in the d = 2 case that ensures -P(i|i) coincides with our normal 
definition of the Pauli Y = cr y operator. The Pauli group Q for a single particle is comprised of the Pauli 
operators along with a set of global phases, 


I {i kp (x\z) | k e Z 4 , x, z e Z 2 } (d = 2) 

\ {u k P( x \z) | k,x,z e Z d ) (d > 2) 

Using so-called symplectic notation, the general form for an n-particle Pauli operator is 


P(x\z) = (* X1 ® X" 2 ...) (Z* 1 ® Z^ 2 ...) (4) 

where x and z are now vectors of length n i.e., x, z e Z d , where each x, or z 3 is an element of Z d . Two 
operators P( x \ z ) an d P( x '\z') commute if and only if the symplectic inner product between their vector 
descriptions vanishes i.e., if and only if Yhi{ x i z i ~ x \ z i) = 0- 

Because of the graph-based contextuality formalism, we are particularly interested in the projectors 
associated with the Pauli operators. The c o k eigenspace of a single-qudit operator ^Wi) corresponds to 
the projector 

nfaubOIfc] = (5) 

is z d 

which is clearly a qudit stabili z er state i.e., it is a rank-1 stabili z er projector. There are a total of d{d + 1) 
distinct single-qudit stabili z er states [10] in a Hilbert space of dimension d, and they can be identified as 
the eigenstates of the following set of operators 


MUB operators: {P (0 |i), P(i|o), P(i|i), • • •, P(i|d-i)} 


( 6 ) 
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Measuring a two-qudit Pauli operator corresponds to projecting with the following rank-d operator 

n(xi,x 2 |2l,Z2)[fc] = y y ^ ^ (-f(a:i,a:2|zi,^2))' 7 (7) 

This projector can be decomposed into a sum of d rank-1 projectors via, 

n■(x 1 ,x 2 |z 1 ,Z2)[fc] = y 1 (n(xi|zi)[a] <8> n ( * 2 | 22)[b ]) (8) 

where there are d solutions (a, 6) for a given k in the summation. To prove this last identity, one must 


insert the definitions given in Equation (5) and Equation (7) and use the fact that 

u -{ma+nb) A m ® B n = dj^ u~ nk (A ® B) n since u a{m ~ n) = d5 m , n (9) 

n aeZ d 

a-\-b=k 

The number of distinct two-qudit stabilizer states (see e.g., [8]) is 

{n} tot = d 2 (d 2 + l)(d + 1) (e.g., d = 2 : 60, d = 3 : 360, <7 = 5: 3900) (10) 

Of these states, the number of separable two-qudit states is 

{n} sep = [d(d+ l)] 2 (e.g., d = 2 : 36, d = 3 : 144, d = 5 : 900) (11) 

which implies that the number of entangled two-qudit states is 

{n}ent = d 3 (d 2 - 1) (e.g., d — 2 : 24, d — 3 : 216, d = 5 : 3000) (12) 


Later, we will be interested in whether these sets of stabilizer states exhibit state-independent 
contextuality. 

2.1.2. . The Clifford Group 

The set of unitary operators that map the qudit Pauli group onto itself under conjugation is called the 
Clifford group, C d , 


C d = {CeU(d)\CP m C 1 kP 1 ^} 

where the proportionality symbol denotes equality up to a phase of the form u> k (or i k in the case of 
qubits). The number of distinct Clifford gates for a single qudit system is \Cd\ = d 3 (d 2 — 1), and this can 
be seen by noting the isomorphism 


Cd — SL(2, Z^) k 


established by Appleby [11]. If we specify the elements of SL(2, Z d ) and Z 2 d as 


F 


a ft 

7 5 


eSL(2,Z d ) 



(13) 


(14) 
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then [11] provides an explicit description of the unitary matrix C(p\u) £ C/ in terms of these elements i.e.. 


C(F\u) P(u\ \u2) Up 

Uf = (jf E'to \j)(k\ p? 0 

lE£=i r “ 7fc2 K)( fc l P = ° 


where r = oj 2 1 if d > 2 and r = (—i) if d — 2. 

The canonical two-qudit Bell state, 


l*> 


1 

\fd 


j=d— 1 


3=0 


(15) 

(16) 


(17) 


is the unique +1 eigenstate of d 2 mutually commuting Pauli operators, powers and products of X 0 X 
and Z 0 Z z.e., 

|$)($| = ^ (/ 0 / + X 0 X + ... + (X 0 xy- 1 ) (I 0 / + Z 0 Z- 1 + ... + (Z 0 Z- 1 )^- 1 ) (18) 
A Jamiolkowski state, | J v ) e C' /J , corresponding to a unitary operation (7 G U(d) is defined by 


|Jtr) = (I<8»C7)|^) 


(19) 


This correspondence between operations and higher-dimensional states is known as the Jamiolkowski 
isomorphism. If the unitary, U , used in Equation (19) is a Clifford operation, then the resulting 
Jamiolkowski state is a bipartite entangled stabili z er (BES) state. In fact, any BES state must be a 
Jamiolkowski isomorph of a Clifford gate so that the number of BES states is 

|{n} ent | = \C d \ = d 3 (d 2 - 1) (20) 


We will use this group structure associated with BES states later, when we discuss whether the set 
| (njent| exhibits contextuality. 


2.2. Quantum Contextuality and Graph Theory 

The following argument, as put forward by Mermin [12], helps illustrate the counterintuitive nature 
of quantum contextuality. Consider a set of mutually commuting operators {A, B, C ,...} and attempt 
to ascribe to each measurement a value {u(A) e A (A), v(B) G A (B), v(C) G A(C),...} respectively, 
where A (A) denotes the spectrum of A and so on. If some functional relation / exists, such that 

f(A, B,C,...) = 0 (21) 

is an operator identity, then we should also have that 

f(u(A),u(B),u(C),...) = 0 (22) 

Explicitly providing sets of operators for which this is not true shows that our assumption of assigning 
zz(A) as the outcome of measurement A etc. was unjustified. In other words, we cannot assume that 
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quantum measurement is a process that just reveals pre-existing properties of a state, independently of 
any other compatible measurements that are being performed. 

Arguably, the canonical example of quantum contextuality is provided by the Peres-Mermin magic 
square [12,13]. Here, the mutually commuting sets of operators are triples of two-qubit Pauli operators. 
There are 6 sets of commuting operators, where each set corresponds to a row or column of Table 1. 
Define the operator comprised of the product of all three operators in Row i or Column j as Ii, and ir¬ 
respectively. Each individual two-qubit Pauli operator has eigenvalue ±1, whereas it can be verified that 
Ri — 1 4 — —Cj for all i,j G {1,2,3}. This forces us to accept that v(Ri) = 1 = —v(Cj). The only 
way v(Ri) = 1 could hold is to have an even number of —1 assignments in each row. The only way 
i '(Cj) = — 1 could hold is to have an odd number of —1 assignments in each column. These last two 
statements are mutually contradictory—no such assignment can be found. 

Table 1. One example of a Peres-Mermin (PM) magic square construction provided by 
Aravind [14]. The operators in each cell of this table are two-qubit Pauli operators. If we try 
to assign ±1 values (i.e., measurement outcomes) to each operator in a consistent way, we 
are forced into a contradiction. 



c x 

c 2 

C 3 

Ri 

X ® Y 

Y®X 

Z ® z 

R'i 

Y®Z 

Z®Y 

x®x 

r 3 

Z®X 

X <Z Z 

Y ® Y 


Cabello, Severini and Winter [15] have recently introduced a graph-theoretical generalization of these 
Kochen-Specker type constructions, inspired by an earlier result of Klyachko et al. [16]. Consider a set 
of binary yes-no tests, which we quantum mechanically represent by a set of rank-one projectors, II, with 
eigenvalues A(II) G {1, 0}. Compatible tests are those whose representative projectors commute, and a 
context is a set of mutually compatible tests. Contradictions can be found if we try to assign outcomes 
to these tests independently of their context. 

Non-locality is a special instance of contextuality wherein compatibility of tests is enforced by spatial 
separation. The two-qubit CHSH inequality provides an example of a non-contextuality inequality, 
which in this case rules out a local non-contextual hidden variable theory. If we drop the requirement of 
spatial separation, then a single qutrit suffices to rule out non-contextual hidden variable theories, as we 
now discuss. 

Take a set of five yes-no questions {n 0 ,IIi... ,n 4 } such that [n^II^i] = 0 where addition is 
performed modulo 5. We further impose, as usual, that commuting rank-1 projectors cannot both take 
on the value +1 i.e., the respective propositions are mutually exclusive and cannot both be answered in 
the affirmative.Graphically, we can represent this scenario as a pentagonal graph T (i.e., a 5-cycle) where 
vertices correspond to tests, and adjacent (connected) vertices correspond to compatible and exclusive 
tests. Define the operator E r to be 

w = E n 

ner 


( 23 ) 
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so that in the classical (non-contextual) case £ r counts the number of yes ( i.e ., +1) answers to the tests. 
The maximum number consistent with the rules outlined above is 

(^r)max V = 2 (24) 

Quantum mechanically we can achieve an orthogonal representation of this graph using five vectors 
|V’i) G C 3 (see Figure 1) so that [n^IL,] = 0 = 0 i.e., n, ; = It was shown 

in [15,16] that the maximum value of Ep achievable by this qutrit representation is equal to the maximal 
achievable by arbitrary-dimensional quantum systems, i.e., 

<E r }£“ = i/5 » 2.236 (25) 

To see this, take the maximizing eigenvector (V’max) °f the Hermitian operator E r and note that 

Tr(£rl(OW'L<l) = (£r)S“ = ^5 (26) 

The operator E r can be interpreted as a witness for contextuality, in the same way that Bell inequalities 
are witnesses for nonlocality i.e., 

Tr(E r p) > 2 =>- p exhibits contextuality (27) 

in the preceding 5-cycle example. 

Figure 1. The KCBS contextuality construction [16] involves five projectors of the form 
IT,; = |'0j)( , 0i| e 'K.i, where the un-normalizcd versions of |Q, ; ) and their mutual orthogonality 
relations are given in the graph depicted (connected vertices correspond to compatible and 
exclusive tests). The state IV’max) is maximally contextual, insofar as it maximally violates 
the noncontextuality inequality (Er) ^ v v < 2. 



In the same way that nonlocality can be generalized to a general non-signalling theory (the most 
famous example being the non-local (Popescu-Rohrlich [17]) box construction), contextuality can be 
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broadened to include generalized probabilistic theories (GPT). Once again, the quantity E r can be 
maximized in a meaningful way e.g., in the preceding 5-cycle example we have 

(Sr)S« = \ (28) 

The remarkable result of [15] is that the quantities (Er)max f° r X G {NCHV,QM,GPT} 
and for a general orthogonality graph correspond to well-known graph-theoretical quantities i.e., 
(ct(r), $(r), a*(T)} respectively. In the 5-cycle example these quantities equate to values given 
in Equations (24), (25) and (28). We can use the many relationships between these and other 
graph-theoretical quantities e.g., 


a(T) < i?(r) < a*(T) (29) 

which is true for all graphs, and this implies 

<Sr>S™ v < (S r >S“ < <Sr>SS (30) 

which makes intuitive sense. Graphs for which a(T) < i?(T) indicate that appropriately chosen 
projectors II can reveal quantum contextuality. It is important to note that not all realizations of T 
using a set of projectors will reach the maximum value $(T). Instead, $(T) is an upper bound that can 
always be reached, by optimizing over the set of projectors used and by optimizing the state p whose 
expectation value is (E r ) QM = Tr(pE r ) with respect to the operator £ r . Sadiq et al. [18] discuss 
this non-optimality in the context of nonlocality, where they argue that spacelike separation is not fully 
captured by the orthogonality graphs as we have described. We will see similar manifestations of this 
non-optimality in Section 3.2 where the expectation value (Er) is significantly less than the most general 
upper bound given by i?(T). 

One significant difference between nonlocality and general contextuality is that arbitrary (even 
maximally mixed) states can exhibit contextuality when the set of tests comprise an example of 
state-independent contextuality (SIC), e.g., the projectors associated with the Peres-Mermin magic 
square [12,13]. A necessary but not sufficient condition [19] for a collection of tests {11} to exhibit 
SIC is that x(r) > D, where (I ) is the chromatic number (explained in Section 2.2.1. ) of the 
orthogonality graph T and D is dimensionality of the Hilbert space projectors i.e., II e Tin- 

2.2.1. . Graph Theoretical Prerequisites 

In order to prove most of our claims we will need to appeal to some well-known facts and results 
from graph theory. Here we provide the necessary background. A graph, T, is a mathematical structure 
consisting of a set of vertices, and a set of edges connecting the vertices. If two vertices g and h are 
adjacent (connected) then we denote that as g ~ h. Since we only ever consider undirected graphs, 
this is equivalent to h ~ g. The graph complement T of a graph T is obtained by replacing edges 
with non-edges and vice-versa. An independent set of a graph is a set of vertices, no two of which are 
adjacent. The independence number a(T) G N is the size of the largest such set. The clique number 
cu(r) e N is the size of largest set of vertices in which every member is connected to every other 
member. The clique number of any graph is equal to the independence number of its complement and 
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vice versa i.e., a;(T) = a (I ). The Lovasz theta number $(T) G M is the solution of a certain semi-definite 
programming problem while the fractional packing number a*(T) G Q is the solution of a certain linear 
program (see [15] for details). The integer x(T) G N—the so-called clique cover number—represents 
the minimum number of cliques needed to cover every vertex of T. It is well-known that this serves as 
an upper bound for the fractional packing number i.e., a* (T) < ;\ (T). The vertex coloring problem for a 
graph involves assigning a color to every vertex in such a way that adjacent vertices cannot be assigned 
the same color. The minimum number of colors required to do this is the chromatic number x(T). A 
little thought reveals that the chromatic number is bound below by the clique number i.e., w(r) < x(r). 

An undirected Cayley graph T(G, T ) with an associated finite group G and set T C G is the graph 
whose vertices are the elements of G and whose set of edges is {g ~ h\g~ l h G T}. For an undirected 
Cayley graph without self-loops we must have / ^ T and T” 1 = T. The resulting graph T(G,T) is 
regular i.e., each vertex is adjacent to |Tj other vertices. 

Given two graphs G and H with respective vertex sets {g, g ',...} and {h. h !,...}, we can define 
various graph products whose vertex set is the Cartesian product of {g, g',...} and {h, h',...} but where 
the condition (g, h ) ~ (g 1 . h') can be defined in a number of different ways. For example the OR product 
was already used in [20] in the context of quantum contextuality 

G OR H : (g,h) ~ (g',h') <*=> g ~ g' oxh~ti (31) 

There are also distinct notions of graph addition, the simplest being the disjoint union whereby the 
total vertex set is the union of vertex sets and the total edge set is the union of edge sets. We denote the 
disjoint union of m copies of the complete graph K n as mK n 

3. Results 

3.1. State-Independent Contextuality Using Stabilizer States 

In this section, we will apply graph-theoretical techniques to sets of stabilizer states to see whether 
they are sufficient to exhibit state-independent contextuality (SIC). We begin by concentrating on the 
set of bipartite separable stabilizer states, (n} sep , and show that they are insufficient for exhibiting SIC. 
We next concentrate on the subset of two-qudit stabilizer states that are maximally entangled, {n} en t- 
As we discuss below, this subset is isomorphic to the Clifford group, hence powerful group-theoretical 
techniques can be used to examine the associated orthogonality graph r ent . Finally, we examine the total 
(separable plus entangled) set of bipartite stabilizer states, {n} tot , and show that it cannot exhibit SIC for 
qudits of small odd prime dimension. Tables 2-4 summarize our results concerning the structure of sets 
of projectors and whether they display contextuality. It seems clear that state-independent contextuality 
is much easier to achieve using qubits as opposed to qudits with odd prime d > 2. At a superficial 
level, one could attribute this phenomenon to structural differences in the Pauli group as outlined in 
Equation (3). Because the multi-qubit Pauli group necessitates a dr- th i.e., 4-th root of unity in its phases, 
the structure is somewhat richer and non-contextual contradictions are easy to find. The simpler structure 
of the multi-qudit Pauli group for prime d > 2 appears to make state-independent contextuality more 
difficult, if not impossible, to find. We shall see in Section 3.1.3. that the local hidden variable model 


Entropy 2013 , xx 


10 


provided by an appropriate discretized Wigner function provides a deeper and more general reason that 
forbids state-independent contextuality using stabilizer projectors in certain dimensions. 

Table 2. The order (number of vertices) of the orthogonality graphs associated with 
two-qudit stabilizer states. 



1 Fsep | 

| Tent | 

1 r tot I 

d = 2 

36 

24 

60 

d = 3 

144 

216 

360 

d = 5 

900 

3000 

3900 

d 

[d(d + l)] 2 

d 3 (d 2 — 1) 

d 2 (d 2 + l)(d + 1) 


Table 3. The independence number a of the orthogonality graphs associated with separable, 
entangled and all two-qudit stabilizer states. The quantities with asterisks are those for which 
we do not have a graph-theoretical proof, but which the more general results referenced in 
Section 3.1.3. support. The contextuality formalism of [15] tells us that maximum value of 
(E r ) in a non-contextual hidden variable theory is given by o (T). The values in bold font 
are those that differ with the corresponding quantum expectation value given in Table 4, and 
therefore the associated projectors exhibit contextuality. 



sep) 

^(r ent) 

a (rtot) 

d = 2 

9 

5 

12 

d = 3 

16 

24 

340 

d = 5 

36 

120 

156 

d > 2 

(d + l) 2 

d(d 2 - 1)* 

(d 2 + i)(d + iy 


Table 4. The achievable values of (Er) in quantum mechanics, for each of the relevant 
sets of stabili z er projectors. This happens to coincide with the clique cover number y(T) of 
the related orthogonality graph. Comparing with the previous Table 3, we see that the only 
situations for which quantum contextuality appears, because a(T) < (Sr) QM , is when we 
include the set of entangled two-qubit states. 



<Sr.,) Q “ = x(r, cp ) 

(Sr» t ) QM = x(r«) 

(Sr„) QM = x(r„,) 

d = 2 

9 

6 

15 

d = 3 

16 

24 

340 

d — 5 

36 

120 

156 

d 

(d + l) 2 

to 

1 

I — 1 

(d 2 + l)(d + l) 
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3.1.1. . Separable Stabilizer States 


Associated with each of the d + 1 MUB operators in Equation (6) is a complete orthonormal basis of 
d states. Furthermore, states from different bases always have the same overlap. To identify individual 
states and bases we introduce the following notation 


{£© = {B\, B\, B \,. ■ ■,, B‘, B \,..., Bf } 


(32) 


and the overlap between any two states is given by 

Tr (B a X') = i(l - 5 a , a 0 + 8 a , a '8 b , v (33) 

Two stabilizer states B% and By are orthogonal if and only if they are distinct members of the 
same basis 


Tr {B a b B a y) = 0 0=^ a = a and b ± b' (34) 

In order to examine the orthogonality graph T s i ng i e of these single-qudit stabilizer states, we associate 
with each state a vertex v(a, b ). Clearly there are |r single | = d(d + 1) vertices with edges 

v(a, b ) ~ v(a', b') <0=> a = a' and b ^ b' (35) 

so that T single is isomorphic to d + 1 copies of the complete graph on d vertices, 

r single ~ [d + 1 )Kd (36) 

Maximum independent sets and maximum cliques are easy to write in this notation; 

Maximum independent set: e.g., (u(a, 1)| a G (1, 2,..., d + 1}} (37) 

Maximum clique: e.g., (u(l, b)\ b e (1, 2,..., d}} (38) 


Next consider what happens if we take tensor products of single-qudit stabilizer states with each 
other. There are [d(d + l)] 2 separable two-qudit stabilizer states of the form 0 By , indexed by 
ai, a 2 G {1, 2,..., d + 1} and b\, b 2 G {1,2,..., d}. The orthogonality of these two-qudit states is 
determined by 


Tr 


B a h \ 0 B 


0.2 

b2 



ci\ = a\ and b\ ^ b\ 
0 or 

a 2 = a' 2 and b 2 ^ b' 2 


(39) 


As before, we can associate with each state 0 By a vertex v({ai, b \}, (a 2 , b 2 }) and the resulting 
orthogonality graph, T sep , has edges 


a i — o/, and b\ ^ b\ 


or 


a 2 = o! 2 and b 2 ^ b 2 


u({ai, 6i}, {a 2 , b 2 }) ~ u({ai, b[}, {a 2 , b' 2 }) 


(40) 
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Note that this is exactly the graph that one would obtain by taking the graph product [see Equation (31)] 
of r s i ng ie with itself 


r = r 

<jpn 


sep 


single 


ORT 


single 


(41) 


With our notation defined as above, we can again describe relevant subsets of vertices i.e.. 


Maximum independent set: e.g., (n({ai, 1}, (a 2 ,1})| oq, a 2 £ (1, 2,..., d + 1}} (42) 

Maximum clique: e.g., (n({l, fq}, {1, & 2 })| h, h e (1, 2,..., d}} (43) 

This implies that a(T sep ) = (d + l) 2 whereas cu(r sep ) = d 2 . Furthermore, we can cover all \d(d + l)] 2 
vertices of r sep using (d + l) 2 distinct, non-overlapping maximum cliques—simply by letting cq and 
a 2 vary over all possible values in (1, 2,..., d + 1} (in our example clique of Equation (43) we have 
set cq = a 2 = 1). This last fact tell us that the clique cover number (see Section 2.2.1. ) is exactly 
x(T sep ) = (d + l) 2 . The following chain of inequalities holds for any graph T, 


a(r) < #(r) < a*(r) < *(r) 


(44) 


Combining this with the preceding discussion, we see that the orthogonality graph for separable states 
has parameters 


«(r sep ) = t?(r sep ) = a*(r sep ) = x(r sep ) = (d + 1) 2 (45) 

which implies that the separable subset on their own are insufficient to manifest contextuality. 

3.1.2. . Entangled Stabilizer States 

Here we will argue, using purely graph-theoretical techniques, that the set of entangled two-qudit 
stabilizer states {n} en t does not exhibit state-independent contextuality, when d is an odd prime. We have 
only shown this to be true for primes up to 11 because one part of our argument still relies on brute force 
computation. What follows is not necessarily the easiest way of proving our claim, but our recognition 
of the orthogonality graph T ent as a Cayley graph is novel, and may be of independent interest. 

A Jamiolkowski state, \ Ju), [see Equation (19)] for which U G Cd must be a maximally entangled, 
bipartite stabili z er state. Since the set of bipartite entangled stabili z er (BES) states is isomorphic to the 
Clifford group, we can hope to find some group-theoretical structure underpinning the orthogonality 
relationship between these BES states. Two BES states are orthogonal if and only if their isomorphic 
Clifford gates are trace orthogonal 


(Jc\Jc>) = 0 « Tr(C'C') = 0 (46) 

Note that, by definition, CCC" must itself be an element of the Clifford group. Define T as the subset of 
Clifford operators having trace zero, 


T:={Ce C d | Tr(C) = 0} 


( 47 ) 
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Two bipartite entangled states \J C ) and | J c >) are adjacent in T ent if and only if the product C r< C' e T. 
Crucially, this is exactly the structure given to us by a Cayley graph introduced in Section 2.2.1. . Note 
that similar Cayley graphs have also been used in the creation of MUBs [5]. 

To exploit this relationship we should first examine the structure of the subset T. We begin by defining 
the Legendre Symbol, £ d , as 


1 if x is a quadratic residue (mod d) 

£ d {x) — ^ — 1 if a: is a quadratic non-residue (mod d) 

0 if x = 0 (mod d) 

where a quadratic residue is an integer in lL d that is of the form x 2 mod d for some x e Z d . Next we 
quote a result from Appleby [ 11 ] that relates the trace of the unitary operator C( F \u) to its constituents F 
and u 


( Case 1: 

\£d{a)\ = 1 

|Tr(C (F | n) )| = <j | 4 ( 7 )| VdS Ulfi 


|Tr(C, FW )| = 


dd Ult od U2l o 
( Case 2: 
\WF)- 2)| = 1 


^ = 0 a ^ 0 ) 

(Tr (F) + 2 ) 
(Tr(F)= 2, 7 ^0) 
(Tr(F) = 2 ,7 = 0) 

P 7 ^ 0 ) 

(Tr(F) / 2) 


\£ d (-/3)\Vd5 u ^-x {1 _ a)ui (Tr (F) = 2 ) 


(48) 


(49) 


By examining Equations (48) and (49), with particular attention to the three relevant cases satisfying 
Tr(F) = 2, we find that there are exactly 

|Tj = d{d - l ) 2 + (d 2 - 1) + [d(d - l )] 2 = [d{d - 1) + 1 }{d 2 - 1) (50) 

such C(f\ u ) with Tr (C(f\u)) = 0. The structure of the traceless subset of Clifford gates is further 
simplified by using a result of Zhu [21] that decomposes the Clifford group into conjugacy classes. 
Using the notation that the class [C(^| u )] with a representative element C( F \ U ) has size | [C^p^)] |, we 
have 


(51) 

Vw 1 e{l,2,...,(d-l)/2} (52) 

y Ul G {1,2,..., (d — l)/2} (53) 

where v is a quadratic non-residue i.e., $ x G such that v = x 2 mod d. By using these conjugacy 
classes to enumerate the number of traceless Clifford gates we find 

^ d(d 2 — 1) = [d(d — 1) + 1 }(d 2 — 1) 


C 


c 


c 


[SS]|[S] 

[it] [“o 1 ] 


[i?] [“o 1 ] 


= d 2 - i 


d(d 2 — 1 ) 
d(d 2 — 1 ) 


\T\ = (d 2 - l) + 2 


2 
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which confirms our previous counting argument. 

The connection set T is closed under conjugation so that T(G, T) is a normal Cayley graph. This 
implies [ 22 ] that 


if a(T)ce(r) = \T\ then x(T) = cu(T) (54) 

Clearly, {TI} ent contains a complete orthonormal basis of cl 2 vectors, which implies that u;(r ent ) = d 2 . 
In our case, the implication of finding o(T) = p(p — 1) is that x(T) = D = d 2 and therefore the 
set {n}e n t does not exhibit state-independent contextuality. As is well known, the problem of finding 
a maximum independent set or the chromatic number of a graph is NP-hard in both cases. However, 
currently available implementations for clique-finding, e.g., Cliquer [23,24], seem to perform much 
better than algorithms for vertex coloring, e.g., Mathematica [25]. For all odd primes p G (3, 5, 7,11} 
that we have checked we can find [23,24] independent sets of size 

a(Fent) = p(p 2 - 1) (55) 

so that the set {n ent } does not satisfy the (Tj > I) criterion [19] for state-independent contextuality. 
Only for the qubit case does {IT ent } exhibit state independent contextuality /.<?., 

d = 2 : a(r ent ) = 5, x(r ent ) = 5 (56) 


Given these results, it seems natural to conjecture that d = 2 is the only prime qudit dimension for 
which the entangled set {n ent } displays state-independent contextuality. It would be nice to find a purely 
graph-theoretical proof that state-independent contextuality is impossible for all odd prime dimensions, 
i.e., without resorting to numerics. 

Associated with each row or column of the Peres-Mermin square in Table 1 is a complete orthonormal 
basis of four entangled stabilizer states. For example, the basis associated with Row 1 can be decomposed 
as the following four projectors 


' 

Row 1: < 


i 

4 

1 

4 

1 

4 

1 


4 


(I ® I + x ®Y + Y ® x + Z ® Z) 

(/ <g> / - X <g> y + Y <g> X - Z <g> Z) 

(/ <g> / + X < 8 > y -Y <g> X - Z <g> Z) 

(/ <g> / - X < 8 > y -Y < 8 > X + Z <g> Z) 


n 2 

n 3 

n 4 


(57) 


and similarly for all other rows and columns, eventually producing a set of 24 projectors (II!,..., II 24 }, 
which provide a proof of the Kochen-Specker theorem. That the Peres-Mermin square can be recast as 
a statement about 24 projectors is well known [13]. It seems to have gone unnoticed, however, that the 
particular choice of Peres-Mermin magic square in Table 1 is completely equivalent to r ent in the way 
we have just described. 


3.1.3. . Generic Impossibility of Qudit State-Independent Contextuality 

With reference to Tables 3 and 4 we see that the condition for contextuality, a(r tot ) < (£ rtot ) QM , 
holds for the two-qubit case, but not for the two-qudit case when d = 3 or d = 5. The proofs are 
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purely computational so we omit any further details. The results obtained so far prompt the question 
whether state-independent contextuality is ever possible using odd prime dimensional qudits. In fact 
this question has already been answered in previous works that dealt with a discrete quasi-probabilistic 
representation of odd-dimensional quantum operators. Implicit in the discrete Wigner function (DWF) 
construction of Gross [8] is a local hidden variable model of a sub-theory of quantum mechanics, which 
was made explicit in a more recent work [9]. This sub-theory is defined by restricting the available 
states and effects ( i.e ., operations corresponding to either measurements or transformations) to those 
that are positively represented within the DWF. Most relevant for our discussion is the fact that this 
sub-theory includes all stabilizer measurements and all (mixtures of) stabilizer states. In particular, 
a set of odd-dimensional stabilizer measurements applied to the completely mixed state always has a 
local hidden variable model. Viewed another way, this tells us that state-independent contextuality is 
impossible whenever qudit (with odd d) stabilizer measurements are used. State-dependent contextuality 
using odd d-dimensional stabilizer projectors is still possible, however, so long as the measured state is 
not positively represented. This is the case in the examples presented in the next section. 

3.2. State-Dependent Contextuality Using Bell Inequalities 

In the previous section we sought cases for which a sum of projectors Sr had an expectation 
value Tr(pE r ) > a(T), regardless of the state p that was used. In this section we will construct 
non-contextuality inequalities that are violated by only some quantum states, i.e., only some states are 
contextual with respect to this inequality. We will convert qudit Bell inequalities into orthogonality 
graphs, creating non-contextuality inequalities that can only be violated by entangled states. Starting 
from the famous two-qubit CHSH inequality [26] e.g., 

(B) < 2 where B = X®X + X®Y + Y®X-Y® Y (58) 

we will show how this corresponds to a non-contextuality inequality in the sense outlined in Section 2.2. 
A similar analysis has already been done numerous times in the literature (see e.g., [15,18]). 
Subsequently, we describe the resulting graphs and non-contextuality inequalities that arise when a qudit 
(for odd prime d ) version of the CHSH inequality is used. This is our novel contribution. 

In order to see the graph theoretic interpretation of Equation (58), first rewrite B as follows; 

B = X <g> X + X <g> Y + Y <g> X - Y <g> Y (59) 

= 2 (11(1,11,0,0)[0] + n ( i,i|,0,i)[0] + n ( i,i|,i,0)[0] + 11(1,11,!,!)[!] — 2I 4 ) (60) 

= 2 ( Y ( n (iio)[a] ®n ( i| 0 ) [6] ) + Y ( n (iio)M ®n ( i|i )[6] ) 

a-\-b=0 a-\-b=0 

+ Y ( n (l|l)M ® n (l|0)[fc]) + Y ( n (l|l)M ® n (l|0[b]) “ 2I 4) 

a+6=0 a-\-b= 1 


( 61 ) 
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where the second equality is a consequence of the definition Equation (7), while the final equality is a 
consequence of Equation (8). Combining (B) < 2 with Equation (61) and denoting the sum of eight 
projectors as S rcHSH returns the following inequality 


(B) < 2 (62) 

=► <2£r c „ s » - 4I i) < 2 (63) 

^ (^Tchsh) — 3 (64) 

where the latter is of the correct form for a non-contextuality inequality. 

Each of the resulting eight projector terms in Equation (61) is a stabilizer projector of rank one [27]. 
By checking the commutation relations between projectors {Eli,..., fl 8 } we get a orthogonality graph 
Fchsh as depicted in Figure 2. The relevant graph-theoretic quantities are a(T CH sH) = 3, 7(T CHSH ) = 
2 + y/2 and cd^Tchsh) = 4. The first of these quantities tells us that the relevant non-contextuality 
inequality is 

(Ei-cXS" < 3 (65) 

which is exactly what our rewriting of the Bell inequality in Equation (64) told us. 

Figure 2. The graph T C hsh that arises from the orthogonality relationship between 
the set of projectors (TIi,..., TT 8 } as constructed in Equation (61). The corresponding 
non-contextuality inequality (E r ) < 3 is completely equivalent to the Bell inequality 

(B) < 2. 



In order to generalize the above example to qudits, we first rewrite the d = 2 CHSH inequality in the 

27T2 

following form, where as usual cu = e rf , 

(B) < 2 where B = ^ kA j ® (A), B 0 = X, A 1 ,B 1 = Y ) (66) 

j,k 6^2 

Ji et al. [28] (with a subsequent reexamination by Liang el al. [29]) have shown that the form of 
the Bell inequality in Equation (66) can be generalized to work for qudits too. They introduced a 
natural modification of the sum comprising the Bell operator, along with a careful choice of measurement 
operators e.g., 


Aj = u jU+1) XZ j , B k = w (2 1)2 (‘ 2+2fc) 17 2 lfc 


( 67 ) 
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For d = 3 i.e., qutrits we have 

(B) < 9 where B = u njk A] <g> B n k (68) 

j,kez 3 

With this information we can rewrite 

b = 3(ii(i ) i|o i o)[o] + n ( i ;1 |o ; i)[ 0 ] + n ( i,i| 0>2 )[i] + n ( i i i| 1 , 0 )[o] + 

+n ( i,i|i,2)[o] + n(ii|2,o)[i] + n(i,i| 2 ,i)[o] + n(i j i| 2j 2)[o] — 3i 9 ) 

Breaking up this expression even further into 27 distinct rank-one projectors, via Equation (8), we get a 
10-regular graph r CHSH of order 27 [30], with a(r CHSH ) = 6 and $(T CH sh) ~ 7.098. The corresponding 
non-contextuality inequality is thus 


<£r CHS XL HV < 6 (69) 

which is easily seen to be completely equivalent to the original inequality (£>} < 9 in Equation (68). 
By looking at the spectrum of the operator £ rcHSH , we see that the maximum achievable using these 
projectors is A max (£r CHSH ) = 6.4115, which falls some way short of the (£r CHSH )£^L — 7.0981 that is 
achievable if we are allowed complete generality in our set of projectors (ITi, n 2 , ..., II27}. 

For d — 5 we get a CHSH inequality of the form 

(B) < 35 where B = u njk A™ <g> (70) 

n&L% 

j,keZ s 

where the Bell operator B eventually produces 125 distinct rank-one projectors leading to a 36-regular 
graph Tchsh of order 125, with parameters a = 12, x) ~ 18.09 and A max ~ 13.09. For d = 7 we get an 
operator of the form 

(B) < 84 where B = u njk A] <g> ££ (71) 

n£Zy 

j,kez 7 

which eventually produces 343 distinct rank-one projectors leading to a 78-regular graph with parameters 
a = 19 and A max ~ 19.4112. We were unable to calculate ^(E C hsh) for d = 7 due to memory constraints. 
All these quantities are summarized in Table 5. 

Observe that for all d 6 (2, 3, 5, 7} that we have discussed, we see some structure in the graphs T C hsh 
obtained from qudit CHSH inequalities. Whereas it is quite easy to see that all graphs will have order 
| Tchsh I = d 3 , we conjecture based on the examples here that 

T C hsh is a (2d -1 )(d - l)-regular graph (72) 

In the qubit case, the graph T C hsh is known to be a member of some well-studied families of graphs 
i.e., Tchsh is a (1,4)-circulant graph and a 4-Mobius ladder [20]. It would be nice to see if T CH sh is 
expressible in a similar fashion for odd qudit dimensions. 
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Table 5. Properties of the orthogonality graphs r CHSH that we have constructed by 
decomposing a two-qudit Bell inequality into stabilizer projectors. Applying (£ r )J^ v = 
a(r) to these graphs gives the maximum achievable value of (Sr) over all local hidden 
variable theories. If we are allowed complete freedom in our choice of projectors then the 
maximum value of (Sr) achievable in quantum mechanics is given by $(r). The maximum 
achievable value of (Sr CHSH ) using stabilizer projectors is denoted by A ma x(Sr CH sH)- We do 
not know of any general solution for these parameters, nor were we able to calculate the 
Lovasz theta number of T C hsh for d = 7, and so these quantities are denoted by question 
marks. The final column concerns subgraphs that can be found within T CHS h, as discussed 
in Section 3.3. 



1 Tchsh | 

a (rcHSH) 

^max(^r CH s H ) 

^(Tchsh) 

c 2 k+l C r C HSH 

d = 2 

8 

3 

3.414 

3.414 

k = 2 

d = 3 

27 

6 

6.412 

7.098 

2 < k < 4 

d — 5 

125 

12 

13.090 

18.090 

2 < k < 6 

d = 7 

343 

19 

19.411 

? 

2 < k < 10 

d 

d 3 

? 

? 

? 

? 


3.3. Induced Subgraphs and Equivalence of Non-Isomorphic Graphs 

In this section we discuss two additional features related to orthogonality graphs that 
display contextuality. 

Firstly, we note a result due to Cabello el al. [31], which says that any graph displaying quantum 
contextuality must contain (as an induced subgraph) an odd cycle C 2 k+\ or its complement C 2 k+\ for 
some integer value of k > 1. This is a necessary, but not sufficient, condition for a given graph T 
to exhibit contextuality. In Table II of [31] the induced C 2 k+ i or G 2 k+i associated with well-known 
Kochen-Specker sets are tabulated. Since the generalized CHSH graphs that we have discussed are new, 
we have performed a similar calculation and summarized the results in the final column of Table 5. All 
relevant subgraphs associated with the qubit CHSH inequality (and its corresponding graph T C hsh) are 
depicted in Figure 3. 

The second aspect of graph-based contextuality that we would like to note is the possibility of 
degeneracy or non-uniqueness in orthogonality graphs. We provide an example of two non-isomorphic 
orthogonality graphs that are both equivalent to the same Bell inequality. This non-unique association 
between inequalities and orthogonality graphs is a general feature, and has already been noted elsewhere 
in the literature, e.g., [18], but we find it instructive to provide an explicit example using the qubit 
CHSH inequality. 
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Figure 3. All graphs that exhibit contextuality contain odd cycles C 2 k +1 or their 
complements C 2 k+ i [31]. For the case of the qubit CHSH graph, the only such subgraph 
that can be found is the pentagon C 5 . There are 8 distinct induced pentagons within the qubit 
CHSH graph, as depicted. 



By a direct calculation, one can show that the sum of six projectors given in Figure 4b produces an 
operator £ r £, HSH = (X 0 X + X 0 Y + Y 0 X — Y 0 Y + 6I 4 )/4, so that rewriting the CHSH inequality 
leads to 


(B) < 2 (73) 

=*• - 6I4) < 2 (74) 

=7 (Sr'„) < 2 (75) 


where the latter is of the correct form for a non-contextuality inequality. The relevant graph-theoretic 
quantities are a(r^ HSH ) = 2, i7(r/ HSH ) = a/5 ps 2.236 and cC(F/ hsh ) = 5/2. The first of these 
quantities tells us that the relevant non-contextuality inequality is 


(Si 


„ ) NCHV < 2 
CHSH' max — 


(76) 


which confirms the equivalence with the original Bell CHSH inequality. In contrast to the original graph 
Fc hsh, the maximal quantum value of i)(T/ hsh ) can never be reached using separable projectors, as this 
would violate Tsirelson’s bound. The structure of this alternate graphical representation of the CHSH 
inequality can be understood as the graph complement of so-called 5-pan graph depicted in Figure 5. We 
have not seen this explicit description elsewhere in the literature but it is implied by the techniques used 
in [18]. 

Perhaps this non-uniqueness can be used to simplify, or better understand the structure of, the 
graphical representation of qudit Bell inequalities that we discussed in Section 3.2. 









Entropy 2013 , xx 


20 


Figure 4. Two non-isomorphic orthogonality graphs, with vertices explicitly labelled, that 
are both equivalent to CHSH Bell inequality (X ® X + X ® Y + Y <g) X — Y 0 Y) <2. 
We showed in Equations (59)—(65) how application of the CSW formalism [15] to 
Figure 4a demonstrated this equivalence. A similar argument for Figure 4b is provided in 
Equations (73)-(76) below. Since our claim is that these graphs display contextuality, they 
must include an odd cycle or its complement; hence we highlight the pentagonal graph C$ 
contained within both graphs, (a) The graph T CHSH comprised of 8 projectors; (b) The graph 
Tchsh comprised of 6 projectors. 




Figure 5. The structure of the alternate CHSH contextuality graph r[ HSH of Figure 4b can 
be recognized as the complement of the 5-pan graph. The family of rt-pan graphs consist of 
an n-cyclc in addition to a single vertex that is connected to the //.-cycle by a single edge. 



4. Summary 

We have applied the graph-based contextuality formalism of Cabello, Severini and Winter [15] 
to sets of two-qudit stabilizer states. Sets of states that arise naturally because of structural or 
group-theoretical significance are seen to display state-independent contextuality only if the particles are 
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qubits, as opposed to qudits with prime dimension d > 2. While our use of graph-theoretical techniques 
was insufficient to prove that this should always hold for particles of prime dimension d, we noted a more 
general result based on the discrete Wigner function that implied that this is indeed the case. It was 
already known that the two-qubit Bell CHSH inequality could be recast as a non-contextuality inequality 
with a related orthogonality graph [15]. Here we have transformed a qudit analogue of Bell CHSH 
inequality into an orthogonality graph and derived the related non-contextuality inequality and graph 
parameters. Along the way we have highlighted various structural properties of the orthogonality graphs 
under consideration. 
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